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Abstract 



We present infrared and Raman measurements of magnetite (Fe304). This 
material is known to undergo a metal-insulator and a structural transition 
(Verwey transition) at Ty = 120 K. At temperatures below Ty, we observe a 
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strong gap-like suppression of the optical conductivity below 1000 cm~ . The 
structural aspect of the Verwey transition demonstrates itself by the appear- 
ance of additional infrared- and Raman-active phonons. The frequencies of 
the infrared-active phonons show no significant singularities at the transition 
whereas their linewidths increase. The frequency and linewidth of the Raman- 
active phonon at 670 cm~^ changes abruptly at the transition. For T < Ty, 
we observe fine structures in the infrared and Raman spectra which may indi- 
cate strong anharmonicity of the system below the transition. Our estimate of 
the effective mass of the carriers above the transition to be m* « 100m, where 
m is the electron mass. Our measurements favor a polaronic mechanism of 
conductivity and underline the importance of the electron-phonon interaction 

in the mechanism of the Verwey transition. 

71.30.-Fh, 71.38.-hi , 78.30.-j, 78.30.Er, 78.40.-q, 78.40.Kc 
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I. INTRODUCTION 



A. The Verwey transition 

Charge ordering in strongly-correlated electron systems has recently attracted a lot of 
attention on account of its possible influence on charge transport and dynamics. In the 
high-temperature superconductors, there is evidence for the formation of charge and spin 
stripes near the commensurate doping x = 1/8 in La2_a;Sra;Cu04 in which superconductivity 
is suppressed.0 The suppression was attributed to the pinning by lattice distortions of the 
dynamic charge and spin stripes formed in this phascS Very recently Raman scattering in 
NaV205 compound revealed considerable charge ordering there at T<35 K.i In the man- 
ganates, charge ordering is observed at certain commensurate doping.! The possibility of the 
coexistence of ferromagnetic and charge-order correlations may be relevant in understand- 
ing the phenomenon of colossal magnetoresistance in manganates. The earliest compound 
known! to manifest charge ordering is magnetite (Fe304). 

Fe304 is a ferrimagnet below 850 K. At room temperature Fe304 has an inverse spinel 
cubic structure with the Oj^ (FdSm) space group.! The unit cell consists of 8 formula units. 
There are two types of Fe positions in this structure. The so-called A positions are charac- 
terized by the tetrahedral oxygen surrounding the Fe ions; the B positions have octahedral 
oxygen coordination. The A positions are occupied only by Fe"*"^ ions, whereas the B posi- 
tions are occupied with equal probability by Fe"^^ and Fe"*"^. 

Magnetite undergoes a first-order transition (the Verwey transition) at Ty = 120 K, 
with changes of crystal structure, latent heat, and a two-order of magnitude decrease of the 
DC conductivity. This transition was discovered by Verwey! in 1939; in his model, the low 
temperature phase has the Fe"*"^ and Fe"*"^ ions at the B-sites ordered, forming layers of ions 
of each charge state. Above Ty there is no long-range order, but short-range order persists. 
According to neutron measurements J3 an orthorhombic space groups (either a 1^1^^ {Pmca) 
or CIjj, (Pmc2i)) may describe the low temperature structure. 
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The Verwey transition has been intensively investigated since its discovery, bur remains 
unexplained. Basically, there are several ideas which partly explain experimental results 
although none is able to give a complete explanation of the transition. This incomplete 
theory is probably not surprising in a system where the electron-phonon interaction, the 
electron-electron interaction, and electronic bandwidth are of nearly the same importance. 
We shall mention some basic ideas which could lead to an understanding of the transition. 

As already mentioned, Verwey and Haaymani originally described this transition as an 
order-disorder transition. The ordered state is insulating, whereas the disordered state is 
metallic due to the motion of the extra electron on the B-positions. Later, Anderson! argued 
that short range order (SRO) is crucial for understanding the transition. He demonstrated 
that certain SRO configurations are energetically favorable (Anderson's condition). In this 
picture, long range order is destroyed above the transition although short-range order, sat- 
isfying Anderson's conditions, persists above the transition. The presence of SRO above the 
transition was confirmed by neutron measurements. However, there is still a question why 
an abrupt increase in conductivity is observed in spite of the persistence of SRO above the 
transition.0 

Cullen and Callen0 presented a picture of the transition in which only electron-electron 
correlations are relevant. It is assumed, that an Fe^^ ion can be represented as an Fe"*"^ ion 
plus an extra electron. In this case, the physics of Fe304 is just the physics of the interaction 
of this extra electron with a crystal of Fe"*"^ ions. The interaction can be described by 
the Coulomb nearest neighbor repulsion and hopping integral for the extra electron. The 
competition between these two values determines whether an insulator or metal is obtained. 

In contrast Mott,0'0 Chakraverty,lli and Yamada0 emphasized the importance of the 
electron-phonon interaction. In the Mott0 picture of the Verwey transition, the charge 
carriers are either polarons or bipolarons. Above the transition, the system can be described 
as a Wigner glass, a system where carriers are localized by the disorder-induced random-field 
potentials. When the temperature is raised, some of the polarons dissociate leading to a 
hopping type of conductivity. 



One should note however that an activated conductivity could also be accounted for by 
the opening of a charge gap. Recent photoemission studie^l of magnetite indicate a charge 
gap of 150 meV below Ty. A recent optical measurement on magnetite^ found that the 
main temperature-dependent changes of optical conductivity occur below 1 eV, connected 
with a gap-like temperature-dependant increase of optical conductivity at 140 meV, followed 
by a broad polaronic peak at 0.6 eV. The gap-like behavior was attributed to the opening of 



a charge gap. The authors of Ref. |Tg suggested that the charge gap, determined as the onset 
of the polaronic band in the spectra below the transition, is consistent with the gap value 



found in photoemission studies. In contrast to the statement of Ref. |1^, we must note that 
photoemission measures the transitions to the Fermi level, whereas infrared spectroscopy 
measures an optically-allowed transition from one electronic band to another across the 
Fermi level. This means that the optical gap of 140 meV would correspond to smaller gap. 



perhaps as small as 70 meV, measured in photoemission. Moreover, according to Ref. |T^, 
the gap does not close above transition but rather decreases by only 50 meV, which is not 
consistent with the analysis of the optical conductivity of Ref. [TB|. 

In this paper we present the results of systematic temperature-dependant infrared and 
Raman measurements of magnetite and discuss both lattice and charge dynamics, as well 
possible mechanism of Verwey transition. 



B. Optical phonons in magnetite 

Above the Verwey transition temperature, T > Ty = 120 K, magnetite has a cubic- 
inverse-spinel structure belonging to the space group OKFdSm). The full unit cell contains 
56 atoms but the smallest Bravais cell contains only 14 atoms. As a result, one should 
expect 42 vibrational modes. Group theory predicts the following modes for this phase. 

A,g + Eg + Ti, + ST^g + 2^2^ + + 5ri„ + 2T2„. 

The Tig, A2U, Eu, and are silent. Thus, there are five Raman- active modes {Aig + Eg + 
3T2g) and five infrared- active modes (5Ti„). 
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When the system undergoes the Verwey transition, the symmetry of the crystal lowers 
to orthorhombic. The unit cell doubles along the z direction and the diagonals of the 
former cubic unit cell become the faces of the low-temperature unit cell. As mentioned in 
Ref. 1^, it has not been possible to get a complete refinement from the neutron scattering 
data. Nevertheless, we can analyze the proposed space groups and compare the results of 
the group theory analysis with our data. Let us start with Pmca which is a nonstandard 
setting of the D}^,Phcm (No. 57). The unit cell contains 56 atoms, and in this case we 
cannot reduce the unit cell as was possible for the room temperature crystal structure. The 
expected number of modes is 168. Group theory predicts the following modes for this phase 

2?>Ag + + 2AB^g + 195i„ + 16^23 + 27^2, + l^B^g + 2653„ 

The Au modes are silent. Thus, there are 78 Raman-active modes {2'iAg + 2ABig + 16-823 + 
IbB^g) and 72 infrared-active modes (19i?iu + 27i?2« + 2QB^u)- It is clear from this analysis 
that one should expect a dramatic increase in the number of modes in the low temperature 
phase. Some of the new modes are the former Brillouin-zone-boundary modes which become 
visible because of the unit cell doubling; others are due to the lifting of degeneracy of former 
T and -E-type modes. Some modes persist over the transition. For instance, the Aig mode 
of the cubic phase becomes an Ag mode of the orthorhombic phase, having the same Raman 
tensor. 

The second proposed symmetry is C|^, Pmc2i, (No. 26). This group does not contain an 
inversion center, leading to very important consequences for the Raman and infrared spectra. 
Namely, if magnetite has such a structure below Ty, both infrared-active and Raman-active 
type modes should be detected by both spectroscopes. Experimentally, this is not the case 
and we can rule this symmetry out. 

The symmetry of the Raman modes can be figured out by their polarization selection 
rules. The Ag and Eg-tyTpe modes have diagonal Raman tensor so that they are seen only 
when the polarizations of the incident and scattered light are parallel. In contrast, the T2g, 
Big, B2g and B^g modes have off-diagonal components so that they are observed when the 
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polarizations of the incident and scattered light are perpendicular. 

II. EXPERIMENTAL 

The magnetite single crystals used in this work were grown by a chemical vapor transport 
technique using stoichiometric Fe304 microcrystalline powder obtained by reduction reaction 
of ferric oxide (Fe203 ). This procedure yielded near-stoichiometric crystals with typical size 
of 4x4x2 mm^. X-ray diffraction confirmed the spinel-type structure of the crystals. Trans- 
port measurements found the abrupt increase in resistance below T = 120K, characteristic 
of the Verwey transition. Optical measurements were carried out on the optically smooth 
surfaces of the as-grown single crystals. 

The reflectance was measured in the frequency range of 50-45000 cm~^ (6 meV-5.5 eV) 
at several temperatures between 40 and 300 K. The sample temperature was maintained by 
mounting the crystal on the cold- finger of a He- flow cryostat. To cover a wide frequency 
range, we used three optical set-ups. A Bruker IFS 113V spectrometer was employed for 
the far infrared and mid infrared regions (50-4500 cm~^; 6-550 meV), a Perkin Elmer 16U 
grating spectrometer for the near infrared and visible region (3800-24000 cm^^; 0.47-3 eV). 
We observed no temperature dependence of the spectra above 20000 cm~^. Therefore, in 
order to extend the data to higher frequencies (20000-45000 cm~^; 2.5-5.5 eV) we merged the 
spectrum at a given temperature with the room temperature reflectance spectrum obtained 
from a Zeiss grating spectrometer coupled with microscope, and finally in the frequency 
range above 45000 cm~^ we merged our data with the data of Ref. |l^. 

In order to analyze the optical properties of the sample we performed Kramers-Kronig 
transformation of reflectance data. For the temperatures above the transition, we use Drude- 
Lorentz model fit in the low frequency part of the spectrum. Then we calculated the re- 
flectance and used the result for the low-frequency extrapolation. Below the transition, we 
assume constant DC conductivity at zero frequency as a low frequency approximation. For 
the high frequency approximation, we used a power law {uj'"^) extrapolation. 
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The polarized Raman scattering experiment was conducted a wide range of temperatures 
(T = 5-300 K) using a triple spectrometer equipped with a liquid-nitrogen cooled CCD 
detector. The 514.5 nm line of an Ar+ laser was used as excitation with nominally 25 mW 
incident on the sample. The spectra were taken in the back-scattering geometry with the 
scattered light polarized vertically while the incident light is polarized to select a particular 
component of the Raman scattering tensor. 

III. RESULTS AND DISCUSSION 

A. Infrared and Raman-active Phonons 

We now proceed to analyze the infrared and Raman data keeping in mind all of the 
above-predicted changes. The infrared reflectance and optical conductivity are shown in 
Fig. 1 and 2, while the Raman-scattering spectra of Fe304 are in Fig. 3. Oxygen phonon 
modes usually have their frequency above 200 cm~^, whereas vibrations of the heavier iron 
atoms should be at much smaller frequency. In Fig. 1 and 2, two Ti^ oxygen modes are 
found near 350 cm~^ and 560 cm~^. These peaks are similar to the infrared- active phonons 
reported earlier by different groups.0^ Four Raman-active phonons are observed at room 
temperature in Fig. 3. The strongest peak at ^670 cm~^ is either the Aig or the Eg mode 
since it is only present in the parallel polarization geometry (XX). Since we expect the Aig 
mode to have the highest frequency as it involves the stretching vibrations of the oxygen 
atoms along the Fe(A)-0 bonds, the 670 cm^^ phonon is assigned to the Aig mode. In Fig. 
3, Raman spectrum in which the incident and scattered light polarizations are perpendicular 
contains the predicted three T2g modes at 193 cm^ , 308 cm^ , and 540 cm" -\ Verblei made 
a similar assignment of the Raman-active phonons in Fe304 except for the 193 cm~^ which 
he did not observe. Instead, he reported that the lowest frequency T2g mode is at 300 cm~^ 
as observed in the spectrum at T = 77 K. However, this feature seen by Verble at T = 77 
K (< Tv) is related to the structural symmetry-breaking that occurs below Ty as discussed 
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below. 

The Verwey transition manifests itself in both infrared and Raman spectra by remarkable 
changes in the phonon spectrum as the temperature is lowered through Ty. In the Raman 
spectrum, one can see at least 17 new modes below Ty. In addition to the new phonon modes 
there is a broad (about 150 cm~^, full width on half maximum) peak with a maximum at 
350 cm~^. Unfortunately, with so many additional phonon modes it is difficult to judge 
whether this broad peak is a superposition of phonon modes or a background on which 
phonon modes are sitting. 

There are several additional infrared modes as well. In Fig. 1 and 2b, in addition to 
the infrared active phonons near 350 cm~^ and 560 cm~^, weak phonon features appear 
below 300 cm^^. Similar effects are observed in the Raman spectrum at T = 5 K as 
illustrated in Fig. 3. Moreover, the Verwey transition also leads to renormalization of the 
phonon frequencies and linewidths. The infrared- and Raman-active phonons which persist 
over the transition were fitted to Lorentzian lineshapes in order to obtain their frequency 
and linewidths. The results of this analysis are shown in Fig. 4. For the infrared active 
phonons at 350 cm~^ and 560 cm~^, the phonon frequencies increase upon decreasing the 
temperature but without visible anomalies around Ty. On the other hand, their linewidths 
manifest an abrupt increase near Ty. Additionally, for the 350 and 560 cm^^ phonons there 
is a reproducible fine structure in the refiectance; see Fig. lb. At temperatures below the 
transition the lineshape of the 350 and 560 cm~^ phonons can not be described by a single 
Lorentzian, which leads to bigger error bars for the linewidth and frequency of these phonons 
below the transition. See Fig. 4. Note that in contrast to the Chakravertylli result, within 
our experimental error we do not observe any soft optical mode. 

The observed phonon anomalies can be explained in the following way. Above the tran- 
sition most of the phonon modes are highly degenerate modes. When magnetite undergoes 
the Verwey transition, the crystal symmetry reduces to orthorhombic, lifting the degeneracy 
of the modes. We suggest that splitting, determined by the degree of distortion associated 
with the symmetry changes, will not be too large. According to the neutron dat J the Fe-0 
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distances for both A and B sites are generally within the estimated standard deviations 
for the prototype structure; therefore, one expects the splitting to be small, appearing as 
a broadening of the modes above the transition. This is what has been observed in our 
experiment for the two infrared active modes at 350 and 560 cm~^. 

The frequency and linewidth of the Raman-active phonon at 670 cm~^ in Fig. 4 shows an 
abrupt increase in frequency with a concomitant precipitous drop in the linewidth at Ty. The 
symmetry of this mode in the cubic phase is Aig] it corresponds to the in-phase vibrations of 
the oxygen ions forming tetrahedra surrounding the Fe ions in the A positionsB At the same 
time, a similar type of displacement for the cubic-orthorhombic structural transition has 
been found from neutron measurements.0 Therefore, it is plausible to assume that 670 cm~^ 
phonon is directly coupled to the structural displacements at the Verwey transition. This is 
the reason why this mode is so sensitive to the transition. 

B. Charge Dynamics 

At room temperature, in addition to the phonon features, one notes a gradual decline 
of reflectance (Fig. la) towards higher frequencies with a distinctive minimum at about 
10000 cm~^ and some weak and almost temperature independent transitions occurring in 
the visible and ultra violet region. 

With decreasing temperature, the low frequency reflectance (below 1000 cm^^) decreases. 
Below Ty a minimum forms at ~ 750 cm~^, followed by increased reflectance, with a 
temperature- dependant hump around 4000 cm~^. Above 20000 cm~^, the reflectance is 
practically temperature independent. 

The optical conductivity derived from the reflectance data is shown in Fig. 2. Aside 
from the phonon peaks, one notes a decrease of the low frequency conductivity (below 
1000 cm~^) with decreased temperature as well as a broad temperature-dependent band at 
5000 cm~^. The intensity of this band increases with decreased temperature followed by 
slight temperature dependent increase of the conductivity towards higher frequencies. The 
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most interesting aspect of the data is the temperature-dependence, which we will discuss 
next. 



According to these calculations of the single-electron density of states, there are 0{2p) 
bands lying about 4 eV below Fermi level which is found in the middle of the Fe{3d) bands. 
Therefore, the interesting physics relevant to the Verwey transition is connected with these 
Fe{3d) bands. Within the Fe{3d) bands one should pay special attention to the two spin- 
polarized bands. Let us recall that there are two magnetic sublattices corresponding to the 
A and B positions in the unit cell, with the spins parallel to each other within one sublattice 
and antiparallel to the spins in another sublattice. From the electronic structure,ii one 
expects a peak in the optical conductivity associated with the optical transition from B site 
to B site. It is clear that such a peak should be sensitive to the Verwey transition because 
the B sites are the sites where ordering takes place. A good candidate for this peak is the 
5000 cm^^ peak, which changes in intensity with temperature. In the framework of the 
Mott mo del0 of the Verwey transition this is the polaronic peak. The increase in oscillator 
strength above 10000 cm~^ is attributed to the oxygen- iron optical transitions, although 
some optically allowed transitions from A to B site may also contribute. 

Another striking feature which should be explained is a gap-like feature in the optical 
conductivity at 1000 cm~^ at temperatures below the Verwey transition. To analyze changes 
of conductivity we estimate the oscillator strength or spectral weight in the low-energy 
region. The effective number of carriers (per Fe"*"^ ion) participating in optical transitions 
at frequency less than u is defined aJ^l 



where m is the free-electron mass, m* is the effective mass of the carriers, Vceii is the unit-cell 
volume, and is the number of Fe"*"^ ions per unit cell. 

One can clearly see in Fig. 5 an increase of A^e// with temperature. For temperatures 
below the transition, we found the maximum change in this quantity at frequencies around 



We will base our discussion on the local density approximation calculations of Ref. |2^. 
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1800 cm~^. In order to compare the temperature dependence of Neff with the DC conductiv- 
ity, we show the A^e//(1800 cm~^), for different temperatures together with the conductivity 
in the insert of Fig. 5. The sharp increase of N^jj closely resembles the very similar increase 
in conductivity at the Verwey transition. A similar result was shown by Park et al.^ 

One should note that there are several ways to interpret the gap-like feature in the 
conductivity. One can assign such a behavior either to an activated mobility of the form 
/ioexp(— ly/z/fc^T), where Wh is the polaron hopping energy, or to the opening of a charge 
gap A. In later case the carrier concentration will behave like noexp(— A/fc^T), leading 
to activated conductivity as well. It is clear that both effects may alter the conductivity. 
However, the agreement of Nf.ff and the conductivity suggests activation of the free carrier 
density. 

Using the definition of N^jj we can try to estimate effective mass of the carriers. Let us 
assume that Ai'e//(1800 cm~^) is purely due to the free carrier optical conductivity. In this 
case, using the definition of A^e//? we can expect A^e//(1800 cm~^, 300K)?a m/m*, which 
gives us m* ^lOOm. Now let us compare this number to the effective mass from Mott 
equation:0 

where hu is the zero-point phonon energy. 

We can get the value of Wh from the polaron peak position. Namely, the peak position 
(~ 5000 cm^^) gives us 2Wb, where Wb is the polaron binding energy. In the simple model 
of an electron jumping between two moleculeslll the hopping energy can be estimated as 
Wh ~ 1/2Wb, yielding Wh = 1250 cm~^. The zero-point phonon energy can be estimated 
to be the frequency of the highest-energy optical phonon, 670 cm^^. These estimates give 
m* ~ 200m, the same order of magnitude that we estimated from N^ff. Note that the Mott 
formula contains an exponent which is very sensitive to both value of Wh and ^hu, therefore 
the fact that we have got the right order of magnitude for the effective mass is consistent 
with the idea of a polaron mechanism of conductivity in magnetite. 
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IV. CONCLUSIONS 



We have presented results for infrared and Raman measurements of magnetite. Both 
Raman and infrared spectroscopic techniques observed phonon anomahes at Ty which are 
attributed to the symmetry and structural changes below Ty. We are able to rule out a 
structure lacking inversion symmetry. 

Our optical data seem to fit to the following description. At the temperatures below 
Ty carriers are localized polarons and we have a strong polaron peak at ^ 5000 cm~^ and 
a gaplike decrease of optical conductivity below 1000 cm~^. All polarons are bound and 
position of the polaron peak occurs at twice the binding energy. 

At temperatures above the transition, some carriers are delocalized, leading to removal 
of oscillator strength from the 5000 cm~^ peak, broadening, and an increase of the low- 
frequency conductivity. It is difficult to separate effects of the filling of the charge gap or 
the activated increase of mobility. However, no matter what is causing it, the low frequency 
optical conductivity follows the DC conductivity. The increase of the 5000 cm~^ band 
oscillator strength with decrease of the temperature is consistent with the polaronic origin 
of this band, consistent with the previous observations of Park et ahS The effective mass 
of the carriers estimated from the optical conductivity gives a value of ~ 100m, close to 
that predicted by Mott in the framework of the polaron model of the Verwey transition. All 
these observations strongly indicate the importance of lattice dynamics effects, favoring the 
polaronic picture of the Verwey transition. 
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FIGURES 

FIG. 1. a) The reflectance of magnetite in the frequency range from 50 to 45000 cm~^ b) 
The low- frequency (up to 700 cm~^) reflectance of magnetite at temperatures above and below the 
Verwey transition. 

FIG. 2. The optical conductivity of magnetite calculated by the Kramers-Kronig transforma- 
tion of the reflectance data. Fig. 2a shows the optical conductivity below 45000 cm~^ and Fig. 2b 
shows the low frequency optical conductivity (up to 700 cm~^) at temperatures above and below 
the Verwey transition. 

FIG. 3. The Raman spectrum of magnetite at temperatures above and below the Verwey transi- 
tion. Shown are spectra in the XX-geometry and XY-geometry. The insert shows the XX-spectrum 
below 600cm~^. 

FIG. 4. The frequency (solid circles) and linewidth (opened circles) of the infrared and Raman 
modes as determined by Lorentzian fits to the data. 

FIG. 5. The effective number of carriers participating in optical transitions at frequencies below 
u {Ngff{u;)) versus u at temperatures above and below the Verwey transition. The insert shows 
Ngff at 1800 cm~^ versus temperature (solid circles), and the DC conductivity versus temperature 
(open circles). 
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